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Abstract 



A sizeable difference in the differential production cross section of top- compared to antitop- 
quark production, denoted charge asymmetry, has been observed at the Tevatron. The ex- 
perimental results seem to exceed the theory predictions based on the Standard Model by 
a significant amount and have triggered a large number of suggestions for "new physics". 
In the present paper the Standard Model predictions for Tevatron and LHC experiments are 
revisited. This includes a reanalysis of electromagnetic as well as weak corrections, leading 
to a shift of the asymmetry by roughly a factor 1 . 1 when compared to the results of the first 
papers on this subject. The impact of cuts on the transverse momentum of the top-antitop 
system is studied. Restricting the tt system to a transverse momentum less than 20 GeV 
leads to an enhancement of the asymmetries by factors between 1.3 and 1.5, indicating the 
importance of an improved understanding of the tt-momentum distribution. Predictions for 
similar measurements at the LHC are presented, demonstrating the sensitivity of the large 
rapidity region both to the Standard Model contribution and effects from "new physics". 
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1 Introduction 



Top quark production at hadron colliders is one of the most active fields of current theoretical and ex- 
perimental studies ITJ. Theoretical predictions fl2l [3l @J [51 for the total production cross section are 
in very good agreement with experimental results both at the Tevatron at 1.96 TeV |[6l 13 and the 
LHC at 7 TeV jH |9). In contrast, sizable differences have been observed between theory predic- 
tions ifTOl [TR [T2l [T3ll for the top quark charge asymmetry and measurements by the CDF and the DO 
collaborations (HI HH HE1 H2l HH 1^3 1211 1221 1^1 at the Tevatron. The discrepancy is particu- 
larly pronounced for the subsample of ti pairs with large invariant mass, m t i > 450 GeV, where a 3.4<r 
effect has been claimed lfi~9ll . It is interesting to note, however, that the discrepancy is less prominent 
in the laboratory frame |fi~9~l . These discrepancies have triggered a large number of theoretical investi- 
gations, using these results, either to restrict new physics like heavy axigluons Il25ll26l or to postulate 
a variety of new phenomena in the t-channel G71 [28l l29l (see also ll30l for a recent review). At the 
same time the robustness of the leading order QCD prediction has been studied in OT1 l32ll . where it 
has been argued that next-to-leading (NLL) as well as next-to-next-to leading (NNLL) logarithmic cor- 
rections do not significantly modify the leading order result, in agreement with the approach advocated 
in lfTTl [i~2l (Note, however, the large corrections observed in Ref. Il33l |34| for the corresponding studies 
of the tt+jet sample). The absence of large corrections in the asymmetry is at variance with the predic- 
tions based on Monte-Carlo simulations where the numerator is evaluated in O(a^) and hence leading 
order (LO), the denominator also in 0{a z s ), corresponding to terms of leading plus next-to-leading order. 
Inclusion of next-to-leading terms in the denominator leads to a reduction by a factor roughly 0.7. A 
small modification of the Standard Model (SM) prediction arises from inclusion of QED corrections. In 
Ref. lfTTl this effect was estimated to lead to an increase of the the asymmetry by a factor 1 .09, in a recent 
analysis ll35l . however, an enhancement factor of 1.2 has been obtained. Obviously this small increase 
of the SM prediction for the asymmetry cannot resolve the discrepancy between theory an experiment 
mentioned above. 

In view of this ongoing discussion, a reanalysis of the SM prediction seems appropriate. In this short 
note we evaluate the QED and weak corrections to the asymmetry, confirming the results of Ref. Il35ll . 
and compare the SM result with the most recent measurements at Tevatron. Subsequently we study the 
effect of a cut on the transverse momentum of the ti system on the asymmetry. A significant increase 
is observed, even for a cut as high as 20 GeV by typically a factor 1.3. This applies both at the parton 
and the hadronic level. Although the implications of this observation for the actual measurement can 
only be studied quantitatively by a (presently not available) full NNLO simulation, this cut dependence, 
nevertheless, may serve as an indication of the sensitivity of the asymmetry on details of the analysis. 

As noted already in lfTTl [i~2l . a charge asymmetry may also be defined and observed at the LHC. Since 
such an effect can only be observed in the small subsample of qq induced events, specific kinematic re- 
gions must be selected where gluon fusion is suppressed and qq- annihilation is enhanced. A particularly 
sensitive observable is the ratio A(Y) = (N(y t > yi) - N(y t < yi)/{N(y t > y { ) + N(y t < yi) with 
fixed average rapidity Y = (y t + yt)/2. For large rapidities, say around 2, the SM prediction for Au(Y) 
amounts up to 0.05 and might well be detected at the LHC. 

The quantity A t t(Y) is also sensitive to physics beyond the SM. We use axigluons as one particularly 
illustrative example and study the sensitivity to amplitudes comparable in size to those suggested by 
recent Tevatron results. 
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2 Amplitudes and partonic cross section 



2.1 QCD asymmetry 

As shown in [fTTl[T2l . the dominant contribution to the charge asymmetry originates from qq annihilation. 
Specifically, it originates from the interference between the Born amplitudes for qq — > QQ (Fig. [TJi) and 
the part of the one-loop correction, which is antisymmetric under the exchange of quark and antiquark 
(Fig. [Tt) (box and crossed box). To compensate the infrared divergences, this virtual correction must be 
combined with the interference between initial and final state radiation (Figs. [Ik Sh). Diagrams with 
triple gluon coupling in both real and virtual corrections give rise to symmetric amplitudes ([TT1 [T2| and 
can be ignored. The corresponding contribution to the rate is conveniently expressed by the absorptive 
contributions (cuts) of the diagrams depicted in Fig |2j 

As a second contribution to the asymmetry we would like to mention "flavor excitation" involving 
again antisymmetric interference terms of different amplitudes, in this case contributing to quark-gluon 
scattering as shown in Fig.|3](amplitudes with triple gluon coupling again don't give rise to antisymmetric 
terms). Flavor excitation hardly contributes to the asymmetry at the Tevatron. At the LHC, however, it 
enhances the asymmetry in suitable chosen kinematical regions, as discussed in Ref. lfT2l . 

Compact analytical results for the asymmetric parts of virtual plus soft radiation (E g < E^ nt ) and of 
hard radiation (E 9 > E^ ut ) can be found in the Appendix of Ref. [1 1 J. 

Let us recall that the color factors corresponding to Fig. [2k and [2b, after averaging over initial and 
summing over final states, are given by 

C a = ^Tr (T°T fe T c ) Tr(T a T c T 6 ) = ^(^L + fL) , 

C b = ±Tr(T a T b T c ) Tr(T*T»T») = ^(^L - fL) • (D 

Without color factors the contributions to the differential cross-section from the two- and three-particle 
cuts in Fig. [2^ and[2t> are related by 

da a (Q,Q) = -da b (Q,Q) . (2) 

The asymmetric piece thus originates from the d 2 abc term and, in leading order, its form is completely 
equivalent to the corresponding QED case. 



2.2 QED asymmetry 

Already at this point we would like to discuss the closely related QED contribution to the asymmetry. 
Let us start with diagram shown in Fig. [4^ following again Ref. ifTTI . The QCD box leads to a color octet 
and color singlet configuration, and the latter interferes with tt production through the photon. A similar 
consideration applies to interference between initial and final state radiation. These two contributions 
are indicated in Fig. [4k by the two cuts, and in combination lead to an additional asymmetric term which 
can be obtained from the QCD asymmetry through the replacement [|TTL \Y2\ 

-> «QED Qt Qq ■ (3) 
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Figure 1: Origin of the QCD charge asymmetry in hadroproduction of heavy quarks: interference of 
final-state (a) with initial-state (b) gluon bremsstrahlung plus interference of the box (c) with the Born 
diagram (d). Crossed diagrams are omitted. 
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Figure 2: Cut diagrams. 



Another QED term originates from the interference between the gluon-7 box with the QCD Born ampli- 
tude. Since gluons and photon are distinct fields, two contributions as depicted in Fig.Hb andHfc arise0. 
Each of these contributes with the factor given in Eq. ©. In total the relative factor between QCD and 
QED asymmetries amounts to 



/•QED _ g a QED Qt Qq _ «QED 36 




. QtQq (4) 

a s 5 



for one quark species. Let us, in a first step, assume identical functional dependence of the PDFs for u 
and d valence quarks in the proton (modulo the obvious factor two) and similarly for antiquarks in the 
antiproton. Assuming, furthermore, dominance of valence quark-antiquark annihilation in tt production, 
the relative contributions of the uu versus dd induced reactions to the cross section have to be weighted 
with the ratio 4: 1 . The QED asymmetry has to be weighted, furthermore, with relative factors /^ ED and 



'These small terms had been neglected in |[T2l . in ifTTI only one of the two had been included. The present result is in 
agreement with |35l 
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(c) 



(d) 



Figure 3: Origin of the QCD charge asymmetry in hadroproduction of heavy quarks through flavor 
excitation. 



/<? ED respectively. The relative QED contribution thus amounts to 



fQED 
J Tevatron 



4/ QED + yQED CKqeD 56 



a s 25 



0.18 



(5) 



at the Tevatron, and thus to an enhancement of nearly twenty percent of the QCD asymmetry, in good 
agreement with the more detailed numerical studies presented below and with the results of 11351 . Com- 
pared to proton- antipro ton collisions the relative importance of uu versus dd annihilation at the LHC is 
shifted from approximately 4 : 1 to 2 : 1, thus reducing / QED to fggg = (2f$ ED + / d QED )/3 w 0.13, 
which is lower than the result of Eq. © by a factor 5/7. The results using standard PDFs are close to 
these values and will be listed in Sect. 13.21 
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Figure 4: Representative diagrams contributing to the QCD-QED interference term. 



2.3 Weak asymmetry 

Weak and electromagnetic interactions are of comparable strength at energies characteristic for the Teva- 
tron and the LHC. Hence, contributions similar to those depicted in Figs.|4ti,|4j5 and|4t with the photon 
replaced by the Z boson should be considered at the same footing. Let us start with the contribution 



4 



depicted in Fig. [5^, where the Z boson is off-shell with virtuality s 3> m|. The result is obtained lfT2ll 
from the photon contribution through the replacement 



QtQq 



{2I t -AQ t s 2 w ){2I q -AQ q s 2 w ) 1 



16 Syy Cyy 



m 2 z /s 



(6) 



with sly and cfy denoting the squares of the sine and cosine of the weak mixing angle, respectively, and 
I q the weak isospin of the relevant quark. Adopting the weighted average similar to Eq. © we find 



n 



weak _ "QED 36 1 



8 „2 
3 b W 



4(1 



8 „2 > 
3 b W , 



3 b W) 



as 5 1 — m 2 z / s 1Q Syy c 



4.4 x 1(T 3 



(7) 



for the contribution Fig. [5^. Note that, as a consequence of the cancellation between up and down 
quarks, and the smallness of the weak coupling, this result is smaller by more than a factor 10 than the 
corresponding photonic result. For proton-proton collision /^ eak is further reduced down to 7 x 10 -4 
and thus is completely negligible. 

It is tempting to estimate the contributions from diagrams [3h and|5fc along the same lines. Indepen- 
dently of any detailed considerations the same compensation between u- and d-quark contributions will 
arise, and in the limit m| <C s, and assuming that final states with real Z radiation are included in the 
sample of tt events, the analogs of Figs. [5}) and[5J; will contribute identically to Fig. [5^, enhancing the 
correction from 0.5% to 1.5%. Alternatively, one may perform an explicit calculation of Figs. \5)p and[5J;, 
allowing for a separation of real and virtual Z boson radiation. 

Let us, finally mention that contributions to the asymmetry involving the squared electroweak ampli- 
tude qq it are of order «q ED /«| and thus at most of 0(1%). Furthermore, if we would include terms 
of order «Q ED /tt| into consideration, terms of order a 2 s «qed m the total cross-section, i.e. electroweak 
corrections to tt production (leading to corrections of comparable size 061 1371 [381 [391 1401 should be 
included as well. Since even the NLO QCD corrections to the asymmetry of 0(a s ) have not been 
evaluated to date, it seems unnecessary to include these C(aq ED ) terms into consideration. 

Indeed, considering the large positive NLO QCD corrections to the cross-section amounting to about 
30% at the Tevatron, it seems plausible to assign a comparable uncertainty to the asymmetry. However, 
in agreement with [[TTl[T2l[3Tl[32l . we shall assume that the central value of the asymmetry will not be 
shifted by the large corrections, in other words, that symmetric and antisymmetric parts of the cross- 
section are shifted by the same factor. 




(a) 
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Figure 5: Representative diagrams contributing to the QCD-weak interference term. 
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2.4 Numerical results at the partonic level 



Results for the angular distribution of the antisymmetric part of the cross section are shown in Fig. [6] (left 
plots) for fixed partonic center of mass energies, and divided by the total qq induced cross-section: 

i(cos£) = ^^^. (8) 

Cqq d COS 9 

QED induced terms are shown separately for u and d quarks, as well as the partonic asymmetry generated 
by qg collisions. The QED generated asymmetries for u and d quarks follows the proportionality factor 
in Eq. ©, amounting to c^qed / a s ^QtQq ~ 0.38 Q q . The piece generated by qg collisions (flavor 
excitation) remains small in the whole kinematic region of interest. 

At this point it is instructive to study the effect of cuts on gluon emission, which lead to an enhance- 
ment of the asymmetry. To understand the origin of this phenomenon, let us in a first step recall that 
the inclusive asymmetry is positive, which can be qualitatively understood from the requirement that 
the final state configuration with minimal modification of the color field is favored, corresponding to 
minimal change of the direction of the color source. On the other hand, if one request a hard gluon in 
the final state, this configuration is more probable for backward going top quarks, thus emitting hard 
bremsstrahlung. Hence, for events with tagged hard gluons a negative charge asymmetry is expected. 
Conversely, one expects a sizable increase of the inclusive asymmetry, if one cuts on events with real 
gluon emission. This is demonstrated in Fig. [6] (right plots) where the dashed and the dotted curves 
are obtained for a cut on the transverse momentum of the top quark pair < P± ax w i m P? ax cho- 
sen differently for different s. We observe an increase of the asymmetry by a significant amount. The 
effect is strongly dependent on the precise value of p™ ax and becomes more pronounced at larger s, a 
consequence of the relatively larger amount of real radiation in the fully inclusive sample. 

The energy dependence of the integrated asymmetry (again normalized relative to the Born cross- 
section a(qq — > tt), 

A = J A(cos6)dcos6 - J A(cos6)dcos§ (9) 

is displayed in Fig. [7J As in Fig. [6l the left plot in Fig. [7J show the pure QCD asymmetry generated 
by qq and qg events, and the QCD-QED mixed asymmetry for u and d initial quarks. A rapid increase 
of the asymmetry is observed in the region very close to threshold, a consequence of the S-wave-P- 
wave interference of the asymmetry. The right plot in Fig. [7J illustrates the effect on the QCD induced 
asymmetry from introducing a cut on the transverse momentum of the top quark pair p". Again a sizable 
increase of A is observed, which depends strongly on the choice of the cut. 

We use the following values for the top quark mass M t = 173.3(1.1) GeV J4TJ, the strong coupling 
a s (M z ) = 0.1184(7) 02, the QED coupling a QED (M z ) = 1/127, and the square of the sine of the 
weak mixing = 0.23. The renormalization scale in Fig. [6] and Fig. [TJhas been set to the partonic 
center of mass energy, fj, = \f§. 
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3 Hadronic collisions 



3.1 Generalities 

The asymmetry can, in principle, be studied in more detail by measuring in addition rriu, the invariant 
mass of the tt system, and its transverse momentum, which is balanced by the transverse momentum of 
an additionally radiated gluon. Controlling both observables together with the asymmetry would provide 
detailed information on the production dynamics. Although such a measurement is only possible with 
large statistics, a first step into this direction has been performed by the CDF collaboration by sepa- 
rating events with invariant mass of the tt-system above and below 450 GeV. However, for a complete 
investigation the statistical precision is far too small. 

In the following we will present in a first step, our predictions for the Tevatron, subsequently for the 
LHC. We will include the update on the electromagnetic and weak corrections, discuss the implications 
of a p*l cut on the asymmetry, and compare our results to the most recent experimental results. In a 
second step we present detailed predictions for the LHC and identify kinematic regions where the charge 
asymmetry could be observed. In view of the indications at the Tevatron for a sizable excess of the 
asymmetry, we then investigate the implications for the charge asymmetry at the LHC in two benchmark 
scenarios beyond the SM. 

Different choices of parton distribution functions of the MSTW2008 set ll43l are used to obtain the 
theoretical predictions of the asymmetry in hadronic collisions, specifically we consider MSTW2008LO 
and MSTW2008NLO. The dependence of the asymmetry on the choice of PDFs is, however, small. The 
factorization and renormalization scales are varied between jx = m t /2 and /i = y/B. This dependence 
gives the bulk of the estimated theoretical error because the asymmetry if proportional to as(/x). The 
variation of the top mass within its experimental error is also considered; its effect on the asymmetry is 
also small. 



3.2 Tevatron 

Let us start with our updated predictions for the Tevatron. Assuming that the rapidities of t and i have 
been measured simultaneously, one defines the asymmetry 

A _ (Y) _ N(y t >yt)-N( yi >y t ) 

At {Y) ~ N(y t > m ) + Wl > Vt) ' (10) 

where Y = (y t + yi)/2 has been fixed. The results as a function of Y are shown in Fig. [8l An almost 
flat A t t(Y) of around 8% is observed. Two versions of the integrated asymmetry have been introduced 
in Refs. lITQl ITTl IT2ll : the forward-backward asymmetry in the laboratory frame 

_ N(y t > 0) - N(y t < 0) _ N(y t > 0) - Nfa > 0) 
lab N(y t > 0) + N{y t < 0) N(y t > 0) + N{y- t > 0) ' 



and the asymmetry in the tt rest frame 



= N{y t > yt) - N(yt > y t ) 
U N(y t > yi ) + N(y t - > y t ) " K } 
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Table 1 : Predicted asymmetries in the laboratory Ai^ and the tt rest-frame A t i at Tevatron. And relative 
amount of uu, dd gg initiated processes. Predictions are given also for samples with the top quark pair 
invariant mass mg above and below 450 GeV. 



laboratory 




m t t < 450 GeV m tl > 450 GeV 


QCD 


0.047 (7) 


0.024 (2) 


0.084 (9) 


QED uu 


0.0094 


0.0047 


0.0174 


QED dd 


-0.0008 


-0.0004 


-0.0012 


weak uu 


0.0011 


0.0006 


0.0021 


weak dd 


-0.0003 


-0.0002 


-0.0005 


SM 


0.056 (7) 


0.029 (2) 


0.102 (9) 


MCFM [19] 


0.038 (6) 


tt rest frame 


Au 


ma < 450 GeV m tl > 450 GeV 


QCD 


0.072 (9) 


0.052 (4) 


0.106(11) 


QED uu 


0.0145 


0.0101 


0.0219 


QED dd 


-0.0012 


-0.0010 


-0.0015 


weak uu 


0.0018 


0.0012 


0.0027 


weak dd 


-0.0005 


-0.0004 


-0.0006 


SM 


0.087 (10) 


0.062 (4) 


0.128(11) 


MCFM [19] 


0.058 (9) 


0.040 (6) 


0.088 (13) 


relative amount 


inclusive 


mu < 450 GeV m a > 450 GeV 


uu 


0.78 


0.76 


0.82 


dd 


0.14 


0.15 


0.11 


99 


0.08 


0.09 


0.07 



Results for both of them are listed in Table[U This Table also lists separately the contributions from pure 
QCD and from QED, the latter again separated for up and down quark. The relative QED contribution 
(QED uu + QED dd) I QCD as obtained from Table [I] is close to 0.18 for all the cases considered and 
thus in agreement with the expectations based on Eq. ©. Similarly, we also give the results for the 
weak terms in the asymmetry, using the approximation m\ <C s. Note that the weak contribution is 
strongly suppressed, even in comparison with the QED piece, and hence can be safely treated in this 
approximation. The relative contribution of gg, uu and dd initiated top quarks is also listed in this Table. 
Due to the deviation of the relative amount of u- and d- contributions from the simple approximation 
4 : 1 we find a slight deviation from Eq. ©. The overall factor 1.21 is consistent with [|35l . 

In view of the strong dependence of the asymmetry on the invariant mass of the tt system, evident 
already from Figs. [6] and [71 we also study the dependence of the integrated asymmetries on a cut on the 
invariant mass of the tt system. Separating the events into two samples with m t f smaller and larger than 
450 GeV respectively, the average value A& = 0.087(10) moves down to 0.062(4) and up to 0.128(11) 
for the two choices. A similar behaviour is observed for the asymmetry Ai^ defined in the laboratory 
frame. 

It is interesting to compare these results with those based on a Monte Carlo prediction lfT9l based 
on MCFM [@4]]. The enhancement factor of the SM result in Table CD compared to MCFM of about 1.5 
is easily understood: a factor 1.2 originates from the inclusion of QED effects, that had been discussed 
in [TT2l and improved in 11351 . Another factor of about 1.3 originates from normalizing with respect to the 
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Born cross-section instead of the NLO result. Since the asymmetric part of the cross-section is presently 
known to LO only we consider the normalization to the LO cross-section more plausible |[TTl[T2ll3Tll32ll . 

For illustration, we compare these theoretical results in the SM with the most recent measurements 
at Tevatron lfi~4l [P71 [i~8l [i~9l as summarized in Table [2l Predictions and results are shown both for A lab 
and A t i and, when available, also split into two samples with m« larger and smaller than 450 GeV, and 
with | Ay | = \y t — yi\ larger and smaller than 1. We note the nearly universal factor ~ 1.5 between our 
result (SM) and the Monte Carlo simulation (MCFM / MC@NLO), which slightly softens the tension 
between theory and experiment. The errors from the choice of PDFs and factorization scale are small, 
the dominant uncertainty arises from varying the renormalization scale and hence the value of as- If 
we would take the difference between LO and NLO prediction for the total production cross-section as 
measure of the theory uncertainty, the error would increase up to ±30%. A graphical illustration of the 
results in terms of the "pull" (measured in standard deviations) is shown in Fig. [9] (errors from theory 
and experiment are combined quadratically. For the results which refer to "reconstruction level" we use 
the MCFM / MC@NLO results, multiplied by a factor 1.5.) The systematic upward shift of all but two 
Tevatron results is evident. The highest discrepancy, as has extensively been discussed in the literature, 
occurs for samples with m t i > 450 GeV and the charge asymmetry defined in the ti rest frame. Also 
shown in this Figure are preliminary results from CMS B51 and ATLAS ll46l with a slight pull in the 
opposite direction. 

Let us now investigate the impact of cuts on hard gluon (and photon) radiation on A t t(Y). The dotted 
and dashed curves in Fig. [8] show the effect of a cut on for values of p™ ax = 10 GeV and 20 GeV, 
respectively. An increase of the asymmetry by more than a factor 1 .5 in the central region is observed for 
the most restrictive choice of 10 GeV, and even a fairly loose p™ ax = 20 GeV modifies the asymmetry 
by up to a factor 1.3. The dependence on Y, the average rapidity of the t and i, is less flat than in the 
inclusive case. Similar enhancement factors are therefore also present in the integrated asymmetries, as 
displayed in Table [3l A fairly similar behaviour is observed if cuts on both m t t and are imposed, 
as shown in the third and fourth column of Table [3j Note that the p± distribution of the tt-system, 
as measured by the DO Collaboration, seems to be at variance with the prediction based on the NLO 
Monte Carlo simulation ||47l (Leading order simulations like PYTHIA Il48l are by construction not able 
to correctly predict the asymmetry.) Considering the fact that in the course of the experimental analysis 
ti events are separated into a sample with four jets (of p± larger than 20 GeV) and a sample with five 
and more jets (with different reconstruction algorithms and efficiencies) it seems important to verify that 
one nevertheless arrives at an unbiased inclusive sample. 



3.3 LHC 

As discussed in flTOl [TO [12]| it is possible to investigate the charge asymmetry also in proton-proton 
collisions at the LHC, exploiting the small ti sample produced in annihilation of valence quarks and 
antiquarks from the sea. As illustrated in Fig. \\0\ production of quarks with larger rapidities will be pre- 
ferred, antitop quarks will be produced more frequently at smaller rapidities. This observation suggests 
to define the cut-dependent asymmetries 

C[VC) ~ N(\y t \< yc ) + N(\ m \< yc ) (13) 
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Table 2: Recent experimental measurements of the asymmetry in the laboratory and in the ti rest 
frame A t i at Tevatron. Results with \Ay\ = \yt — Vt\ larger or smaller than 1 are also summarized. 
Numbers with * refer to "reconstruction level" |fT9l[l4l . the others to parton level. 



laboratory 




ma < 450 GeV 


ma > 450 GeV 


SM (this work) 


0.056 (7) 


0.029 (2) 


0.084 (9) 


CDF EH 5.3 ftr 1 (1+jet) 


0.150 ±0.050 ±0.024 


0.059 (34)* 


0.103 (49)* 


MCFM/MC@NLO |19| 


0.038 (6) 


-0.008 (5) 


0.022 (7) 


ti rest frame 


A« 


m tt - < 450 GeV 


mtf > 450 GeV 


SM (this work) 


0.087 (10) 


0.062 (4) 


0.128 (11) 


CDF ffUl 5.1 ftr 1 (dilep) 


0.42 ±0.15 ±0.05 






CDF HH 5.3 ft" 1 (1+jet) 


0.158 ±0.072 ±0.017 


-0.116 ±0.146 ±0.047 


0.475 ±0.101 ±0.049 


CDF H3 (combined) 


0.20 ± 0.07 ± 0.02 






MCFM H19II 


0.058 (9) 


0.040 (6) 


0.088 (13) 


JJU ||14|| 0.4 ID (i+jetj 


n. 1 qa /'A^ 
u.iyo (oj ; 


U.U/o (4oJ 


ft 1 K /AA\* 

U.l 1 J (OU,) 


MV^^S'rNLU II 1411 


U.IDU (1UJ 


A A1 1 //IX* 
U.U1J (Oj 


A A/1 1 /I 1U 


ti rest frame \Ay\ < 1 


A t t 


m tf < 450 GeV 


m ff > 450 GeV 


SM (this work) 


0.057 (4) 


0.053 (4) 


0.069 (5) 


CDF HU 5.3 ft" 1 (1+jet) 


0.026 ±0.104 ±0.056 






K r t — i"i — in M it 1 Ah 

MCFM H91 


0.039 (6) 






DO El 5.4 ft" 1 (1+jet) 


0.061 (41)* 






MC@NLO 031 


0.014(6)* 






£t rest frame Ay > 1 


Au 


ma < 450 GeV 


ma > 450 GeV 


SM (this work) 


0.193(15) 


0.149 (8) 


0.209 (15) 


CDF E3 5.3 ft" 1 (1+jet) 


0.611 ±0.210 ±0.147 






MCFM nn 


0.123(18) 






DO dm 5.4 ft" 1 (1+jet) 


0.213 (97)* 






MC@NLO US 


0.063 (16)* 







and 



N(\y t \ >y c )-N{\yt\ > y c ) 



A o(yo)~ N(\ Vt \> V c) + NW>vc)' (U) 

which serve to characterize the depletion of top quarks in the central region (A^(yc) > A ( Q t {yc) for 
y c < 0.7 approximately ifTOl ). and their enhancement at larger rapidities (A^(y c ) < A ( j t {y c ) for 
y c > 0.7). Note that we have defined A^(y c ) in Eq. ([□]> with the opposite sign than in Ref. ifTOl 141)1 l50l 
such that both A^(y c ) and A ( ^ tt (y c ) are positive in the SM. The dependence of A™, and A^ on y c is 
shown in Fig.QT]for \/s = 7 TeV (left plot) and 14 TeV (right plot). As one can observe in these Figures, 
A'q 1 is much larger than A™ at large values of the rapidity cut yc- This is because the central region is 
dominated by gluon fusion processes, while the sample with large rapidities has a larger relative content 
of qq initiated events. The statistical significance of both observables is, however, very similar [f5TTl 
because the larger size of the asymmetry A^ with respect to Aq is compensated by the lower rate of 
events at larger rapidities. We consider also the cut-independent charge asymmetries 

= N(A V > 0) - N(A V < 0) 
c N(A V > 0) + N(A V < 0) K 
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Table 3: SM asymmetries in the laboratory A\ a ^ and the tt rest-frame A t i for different cuts in pf. 



laboratory 


Aiab ma < 450 GeV mg > 450 GeV 


p*t < 10 GeV 
pf < 20 GeV 


0.090(12) 0.047(3) 0.161 (16) 
0.076(10) 0.040(3) 0.137(13) 


tt rest frame 


Au m t t < 450 GeV m tI > 450 GeV 


< 10 GeV 
pf < 20 GeV 


0.136(16) 0.097(8) 0.201 (19) 
0.115(13) 0.082(7) 0.171 (16) 



Table 4: SM cut-independent asymmetries A v and A y at different LHC energies. 





n c 




LHC 7 TeV 


0.0136 (8) 


0.0115(6) 


LHC 8 TeV 


0.0122 (7) 


0.0102 (5) 


LHC 10 TeV 


0.0101 (6) 


0.0082 (4) 


LHC 12 TeV 


0.0087 (5) 


0.0068 (3) 


LHC 14 TeV 


0.0077 (4) 


0.0059 (3) 


LHC 7 TeV CMS 051 1.09 fb" 1 
LHC 7 TeV ATLAS BH 0.7 fb" 1 


-0.016 ± 0.030 t^oio 


-0.013 ± 0.026 t$sm 
-0.024 ±0.016 ±0.023 



and 

= N(A y > 0) - N(A y < Q) 
c N(A y > 0) ± N(A y < 0) ' 1 J 



where A^ = \r)t\ — \rft\ and = — |2/f|, which have been used in the recent CMS 11451 andALTAS 11461 
analysis. The SM predictions for the integrated asymmetries are listed Table 0] for different center-of- 
mass energies of the LHC, together with the experimental results for y/s = 7 TeV. Both experiments 
obtain negative asymmetries, although compatible with the SM prediction within uncertainties. 

Top quark production in proton-proton collisions is dominated by gluon fusion, which, in turn, is 
dominant in the central region. Conversly, quark-antiquark annihilation will be more enriched for events 
with tt at larger rapidities (and larger m t i). This suggest to employ the definition of Eq. (fl"0l ), which is 
essentially the asymmetry in the tt rest frame, also for the present case, and concentrate on tt events at 
large rapidities. The prediction for Au(Y) is shown in Fig. [12] for y/s = 7 TeV (left plot) and 14 TeV 
(right plot). Note that Au(Y) is now, by construction, an antisymmetric function of Y. Since most of the 
charge asymmetry is concentrated at large rapidities the statistical significance of any measurement will 
be enhanced, if the sample is restricted to larger rapidities. Let us therefore define the quantity 

A cut lyr \ 

N{y t >yi)-N{y- t >y t ) 
tl { CUt) ~ N(y t >yt)+N( W >y t )> 

where (y t + y£)/2 > Y cut 0. The prediction for A^ t (Y ai t) is shown in Fig. [13] for y/s = 7 TeV (left 
upper plot) and 14 TeV (right upper plot). Note that the prediction includes QED and weak corrections, 
which amount to roughly a factor 1.1. To estimate the statistical sensitivity of such a measurement, the 



■fThe asymmetry (F C ut) of the sample with (y t + J/f)/2 < — Y cut is of the same magnitude but of opposite sign. 
Both samples can be combined to enhance the statistical significance of the measurement by introducing a flip of sign in the 
definition of A c t f (Y cut ) for events with (y t + y { )/2 < -Y cut . 
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cross section for events with tt at large rapidities \ (y t + yi)/2\ > Y cut is shown in Fig. [13] (lower plots) 
as a function of Y cut . Using y/s = 7 and an integrated luminosity of 20 fb 1 as example, more than 
4 x 10 5 events are expected for Y cut = 1. Even allowing for a significant reduction of the sample by 
small efficiencies, the asymmetry A^ t (Y cat ) of more than 2% could be observed by experiment. 

As discussed before, cuts on and m t i may lead to a significant modification of the asymmetry. 
The impact of a cut p*l < 20 GeV on A t i (Y) and A^f 1 (Y cut ) is shown Figs. [L2l and [T3l (upper plots) as 
dashed curves. Numerical results for the integrated pair charge asymmetry with F cut = 0.7 are listed in 
Table[5jfor the inclusive sample, A™ t (Y cut = 0.7), and for subsamples with m t i larger and smaller than 
450 GeV. By definition A^ t (Y cat = 0) = A y c , for which results are listed in Table[4] 

Table 5: SM integrated pair charge asymmetry at different LHC energies for Y cut = 0.7, for the inclusive 
sample, A^f t (Y CVLt = 0.7), and for subsamples with m t i larger and smaller than 450 GeV. 





A^(Y cut = 0.7) m t f < 450 GeV m& > 450 GeV 


LHC 7 TeV 
LHC 8 TeV 
LHC 10 TeV 
LHC 12 TeV 
LHC 14 TeV 


0.0203 (8) 0.0148 (5) 0.0263 (8) 
0.0178(6) 0.0128(4) 0.0224(7) 
0.0142(5) 0.0104(4) 0.0174(5) 
0.01 17 (4) 0.0085 (3) 0.0143 (4) 
0.0100(4) 0.0075(3) 0.0121 (4) 



4 Charge asymmetry beyond the SM at the LHC 

As noted in ifTOl . the asymmetry induced by a "conventional" axigluon G, i.e. with identical axial- 
vector coupling g s for all quarks and assuming m G > 2m t , is negative. This has lead to stringent 
bounds on m G |[25ll . On the other hand the apparent positive excess as observed at the Tevatron has 
lead to numerous suggestions for physics beyond the SM which, however, are difficult to reconcile with 
other experimental facts (For a recent discussion see e.g. 113010 . On the contrary, both measurements at 
CMS [45] and ATLAS fl6l point towards negative asymmetries, although still with large uncertainties 
and compatible with the SM. Independently of these theoretical and experimental considerations it will 
be interesting to investigate the same phenomena at the LHC. As discussed previously in iTTTTl (see e.g. 
Figs. 1 1 and 12) it is possible to identify kinematic regions where qq annihilation into tt is comparable 
or even larger than gluon fusion, and the charge asymmetry is suitable to probe new physics beyond the 
SM in that kinematical regions. We do not try in this section to focus on a particular model that would fit 
better the Tevatron anomaly and then extrapolate that model to the LHC. Rather, we try to investigate the 
power of the charge asymmetry defined in Eq. (fTTI) to discriminate among two of the simplest models 
giving rise to a BSM charge asymmetry, one of them giving a positive excess, and another one leading 
to a negative contribution. 

The most recent measurements at the LHC of the dijet cross-section |[52l [531 impose stringent con- 
straints on axigluon masses below 3 TeV. Still, those limits can be relaxed when considering top quark 
pair production in models in which the coupling of the axigluon vector boson to light quarks is much 
smaller than the coupling to the top quark. As the simplest case, however, we consider here the case 
were the extra gauge boson couples with the same strength to light and top quarks, in two different sce- 
narios: a flavour universal case (octet U), i.e. vector-axial couplings to light and top quarks are equal 
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to the strong coupling g s multiplied by a factor 1.8, g\ = g\ = 1-8, and a flavour non-universal case 
(octet A) with vector-axial couplings to light and top quarks of the opposite sign ll25l . g q A = —g A = 1.8. 
The latter naturally produces positive contributions to the charge asymmetry, and has been advocated as 
one of the possible solutions to the Tevatron anomaly. The former gives a negative contribution to the 
charge asymmetry, and is disfavored by most of the Tevatron measurements, but still compatible with 
some of the measurements within 2a (see Fig. |9j). In both cases we consider an axigluon mass of 3 TeV 
as benchmark model. 

The pair charge asymmetries A t i(Y) and A t t(Y cut ) for both benchmark models and for y/s = 7 TeV 
and 14 TeV are shown in Fig. [14] in comparison with the SM prediction. We also provide predictions 
for samples with a large invariant mass of the top quark pair m t i > 450 GeV As expected, the BSM 
contribution to the asymmetry in the octet U model is negative, and it is positive in the octet A model. For 
low values of F cut the integrated asymmetry A t t(Y cut ) is almost twice the asymmetry in the SM for octet 
A, and almost vanishes for octet U. Introducing a cut in m t i also enhances the size of the asymmetry, as 
in the SM, and in particular for large values of Y cut . 



5 Summary and Conclusions 

The Standard Model predictions for the top quark charge asymmetry have been reanalysed including 
QED and weak corrections corrections. For proton- antipro ton collisions QED terms lead to an en- 
hancement by a factor of about 1.2 in agreement with [|40l . and slightly larger than the factor 1.1 ob- 
tained in [11 J. In total, our prediction is larger than the NLO Monte Carlo results by a factor around 
1.5 ~ 1.2 x 1.3, where the second factor arises from the different normalisation prescription. 

The effect of a cut on the tt transverse momentum has been studied and shown to lead to a significant 
enhancement of the asymmetry. As a characteristic example we study a pj-cut of 20 GeV which leads 
— for the Tevatron — to an enhancement of around 1.3, and even more for more restrictive cuts. 

Various definitions of observables are presented which are sensitive to the charge asymmetry and 
which can be measured at the LHC. The quantity A t t(Y), which measures the forward-backward asym- 
metry with respect to the average rapidity of top and antitop quark, can amount up to 7%, if the large 
F-region is selected. Considering the large statistics expected for the LHC in the near future, this asym- 
metry (and its integrated version) might soon be measurable at LHC experiments. We have also provided 
predictions in two benchmark axigluon-like models for these new observables. 
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Figure 6: Left plots: differential partonic asymmetries from qq and qg induced events for different 
choices of the partonic center of mass energy. The asymmetry from mixed QED-QCD contributions is 
shown separately for u and d quarks. Right plots: differential partonic asymmetries after introducing a 
cut in the transverse momentum of the tt pair p 1 ^. 
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Figure 7: Left plot: integrated partonic asymmetry from qq and qg induced events for different choices 
of the partonic center of mass energy. The asymmetry from mixed QED-QCD contributions is shown 
separately for u and d quarks. Right plot: integrated partonic asymmetry after introducing a cut in the 
transverse momentum of the tt pair pf. 
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Figure 8: Rapidity dependence of the asymmetry A t i (Y) as defined in Eq. (flOl) . Solid line: no cut on p l [, 
dotted/dashed line p^ ax = 10 GeV / 20 GeV. 
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Figure 9: Summary of experimental measurements of the charge asymmetry in comparison with the SM 
theoretical predictions. The histogram represents the pull of the discrepancy for each measurement. The 
asterisk * indicate "reconstruction level" results, the others are parton level. 
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Figure 10: Preferred and suppressed configurations at the LHC. 
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Figure 1 1 : In and out charge asymmetries at the LHC with yfs = 7 TeV and 14 TeV as a function of 
the cut yc in rapidity. The out charge asymmetry A c Jt (yc) is also calculated with an upper cut in the 
rapidity of the top and antitop quarks \y t \, \yi\ < 2.5. 
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Figure 13: Integrated pair charge asymmetry and integrated cross-section as a function of Y cut at the 
LHC with y/s = 7 TeV and 14 TeV. Solid line: without extra cuts, dashed line: p™ ax = 20 GeV, dotted 
line: m t i > 450 GeV. 
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Figure 14: Pair charge asymmetry A t i at the LHC as a function of the mean rapidity Y — (y t + yi)!^, 
and integrated pair charge asymmetry A^f* as a function of Y cut , with and without a cut in the invariant 
mass of the top quark pair m t i. Left plots with y/s = 7 TeV and right plots with 14 TeV. 
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